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ABSTRACT 


In a wedge shaped ocean, the method of 


rmages 1S Usea to 


develop an analytical approximation of the acoustic pressure 


field. Contemporary work develops acoustic 


doublets froma 


combination of the source and surface reflection image using 


simple dipole theory. The method of images 
sum the dipole images. This thesis matches 
achieve a gquadruplet expansion. A computer 
derived quadruplet equation is then created 


results by comparing them with the "URTEXT" 


ere 


1s then used to 
dipole pairs to 
program using the 
to verity the 


program. 
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ls INTRODUCTION 

This thesis 1S a continuation in the examination or the 
acoustic pressure field 1n a wedge shaped ocean. As the 
focus of naval operations shifts to the littoral regions, 
more attention 1S being paid to the problem ot A.S.W. in 
acoustically shallow water, 1.e. regions where the sound 
paths have multiple interactions with the surtace ana bottom. 
To date, the simplest propagation model has been the method 
of images. 

In its present form, this model asSumes 1So0Speed water, 
a pressure release upper surface, and a lossy, penetrable 
Bottom. For parallel surtaces, the source and 1tS images 
form a vertical array with the actual source at the center. 
Furthermore, it tne source 1S near the surface 1t can be 
combined with 1ts negative (igO" cut ot phase) surtrace 
reflection image to form a dipole source. The column ort 
images can then be considered a column or doublet images. 

The equations used in this thesis are the far friela 
approximations of the acoustic doublet. An unbalanced 
doublet is defined as having two sources of nearly equal 
amplitude but opposite phase separated by a distance da. It 
the amplitudes differ greatly, or if the phase difference is 
not 180°; then the doublet approximation 1s not valid. 

When the method of images 1S used 1n a wedge shaped 
ocean, the vertical column ot doublet images becomes a circle 


v 


centered on the apex. In previous work the upper and lower 

image doublets were summed individually with great success. 
This thesis attempts to group upper and lower doublets into 
double doublets or quadruplets. 

This method can only be used tor extremely small wedge 
angles. With larger angles the sound paths from the upper 
and lower doublet images will be to dissimilar, and the 
quadruplet approximation will be invalid. Since typical 
slopes for most of the worlds continental shelves are about 
3°, a good analogy would be that ot using a ship mountea 
active sonar searching for a diesel submarine i1n the shelt 


region. 
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II. WAVE NUMBERS AND SCALING 
Sound propagation ina shallow channel with a pressure 
reiease surface parailei to a rigid bottom has a propagation 


wave number for each mode defined by the following equation 


v15 
H F(6) 7% 


Where H is the depth of the channel. For a fast bottom, 8 


is the critical angle found from 


c 
sind.= = (20 
2 


An analogous reference angle 0. tor a slow bottom 1s 


tan8,-— (3) 


If the depth H 1s replaced by the scaling distance, 


[2D = 


& 


(4) 


& | by 


and the bottom 1S tilted at angle ®; then eq.1 for the 


cutotf value ot the lowest mode with a fast bottom becomes 


obs v1 
2sin@8tanB (4 


and for a slow bottom 


v1 


Kx= 2232S 3 
2tan@tanB ee 


III. ACOUSTIC DOUBLET FORMATION 


The sound propagation for spherical spreading of a point 


eeurce is 


I 


cp) 
When two point sources of opposite phase are combined to form 


a doublet, equation 7 can be expanded, 


Aa ions | 
the oe eslve oe eslor kr,) (8) 


a 


Where the "+" and "-" subscripts reter to the upper and lower 
sources respectively. These values can be expanded 
bE e/a A,=A,+4A 


uy) 
EN 109 A_=A,~AA 


Where the delta values indicate the incremental ditrerence 


tor each source compared to a tneoretical point source 


located exactly between them. The pressure equation can then 


be rearranged to yield 


1+AA 4-44 
Po A2 ei(wt-ks) ae ) @kbr_(__ Ae) @ jhary (10) 
- oa 7+ 22 
i Zz 


An acoustic source located near the surtace has a 


Surface reflection image of equal amplitude anda opposite 


phase. In the far field the image pair is treated as a 


Single doublet source. For a baianced source, the pressure 


equation can be approximated by 


pe Zeivt [ec Ik (r-Ar) ~ G-jk(rehr) } ( De) 


As seen in fig. 2, Ar is 


Gee 
Ar=—sin@ (12) 


OF 


A7y=-r7,Sinysing 
(13) 


combining eqs. 11 and 13, yields the pressure 
p=293 Asin (kr,sinysin@) ef (@e-kr) (14) 
1g 


To develop an equation tor the unbalanced doubiet, the 
amplitude relationships expressed in eq. 9 must be rewritten 


as 








A,ta A, a 
as AA (15) 
So that for an unbalanced doublet 
P=2 jZ5 [ : k . e = iA A k . @ 7 (we-kr) E 
52 ae sin (kr,sin®) -7 a 6os | Sin) re. ( Lies 


(2) 


The distance from the center ot the doublet to the 


receiver can be calculated using the !aw of cosines. 


or 


r*=r;+r;-2r,r,cos (2np+8) wi) 
Where the ranges r,) and r, are measured from the apex to the 
source, and from the apex to the receiver respectively. The 
angle 6 is measured from the surface down to the recelver. 
The range from each successive doublet in a wedge of angle B 


is 


Y,zV¥litr; -27r,7r,cos (2nB+8) (18) 


Where the angle in the cosine term refers to the angle or 
each reflected image; those with a "+" are the upper images, 
oe oe 


and those with a "-" are the lower images. This equation can 


be rearranged to yield 


cy eg 
—+—~ (1-cos (2nf+5) (19) 


J geal eo 5 Nd a 
Gre! ae 


The first order Taylor's series expansion can tnen be used to 


Simplity eq. 19 


rae (X,-r,) 4+—*"2. [1-cos (2nB38)}) (20) 
tn ts 


FIG. 2 
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Lc. QUADRUPLET EXPANSION 
The upper and lower image doublets can be combined to trorm 
a quadruplet. The complimentary doublets are 180° out or phase, 


and the pressure of the new quadruplet 1s 


“A - 
p,--2*(sin [kr,ysin(2nB+8) Je jkAr, 


in [i ¥ 7 ade von 
-gin[kr,ysin (2nB-8)]e2**" Jester 


In this case, both of the interior exponents involving Ar are 
negative. The Ar terms are derived trom a comparison to tne 
range trom the primary doubiet, thereby ensuring that ali pnase 
angles are calculated with respect to a common reterence point. 
To find the Ar values used in eq. 21, an approximation or eq. iv 


gives 


- - Jha Je 
Ari-r!-(z,-r,) = 2° [1-cos (2nB18)) (22) 
‘ci 


The cosine term can be expanded and approximatea tor small 
weoreues Of Oo. 
cos (2nB+4) =cos (2n$) +6sin(2nB) (23) 


Equation 21 can now be expanded into the unbalanced guadruplet 


equation. 





r4z. 
2zA -~jk a 2 {i-cos (2nB) 
2 D w C-kI) a es 
rere ie e one 
ee oo ¢ ‘ 
Cierer écin/2r8) (24) 


{sin [kr ,ysin(2nB+8)}¢ 





, ik =a" 8sin(2nB) 
Zone Vouniznp-0))e “2 "* 


It 1S convenient to define 


b=kzr,p (1-cos (2nB) ) 
d=kr,dsin(2nB) 
Z2 
I,~I, 


(25) 





p= 


So that eq. 24 becomes 


= 2J ej \we-kr) 9-56 
ee e 
(a, (sin [kr,ysin (2nB+6)] -i—cos [kr,ysin (2nb+6) |] e J# (26) 


-B,{sin[kr,ysin(2nB-8)] -j42cos [kr,ysin(2nB-6)]})eae 


The vaiues tor AA and AB must be determined trom tne reriection 


coefficients. 


LOU 


VV. REFLECTION COEFFICIENTS FOR A SLOW BOTTOM 
The reflection coefficients used in this thesis are 


derived from the Rayleigh reflection coefficient 


_ be-sin6,/sin®, 


R= ee a 2s 
be+sin6,./siné, oe” 
Where 
p- Pa c- 2 (28) 
Pp. Cy 


The angles are the grazing angle 68, and angle or transmission 


8, . Therefore, 6. 1s detined as 


(25) 


Using trigonometric relationships, eq. z/ can pe 


rewritten. 


be-/1-c*cos76,/sin6, 
pee = eee (3U) 


be+/1-c*cos*8;/sin®, 
In addition, the term under the radical can be rewritten. 


y1-c*cos78,=c/sin*6,-sin76_ (31) 


Combining eqs. 30 and 31, and rearranging ylelas 


p- Dryi-[sinb./sin8,]” _b-Y1-17x* (32) 
b+y1-[sin8,/sin6,J? b+yi-1/x” 


Weak 


Where x 1S 





sin8, 
aT (33) 
sin8. 
For a slow bottom with a very smail angle, eq. 3z can De 
approximated by 
2be @, 
Ree Vie" a 
The argument of the exponential 1s 
g-_2be .__ 2P2/Pi ee 
yiaC” VC ey) a 
So that this approximation ofr the Rayleigh refiection 
coetficient can be written in amore concise torm 
R=e *: (50) 
Since each ray intersects the bottom and surtrace multiple 
times, a product ot retlection coetticients aiong eacn 
individual path is required. This cumulative coetticient 
will take the rorm 
ie at 
R= []  x{6,-(2n-m) B) (au 
m=1,3,... 
and the product yields 
_, mad; 72 
R=e Wie aos (38) 


For very small angles, this can be approximated py 


i2 


R,=e 2 4h (39) 


The angles of incidence are needed to calculate tne proauct. 


The angles for the bottom images are 
6-2 (u- 1) B (4U) 
Using the law of sines, the angle at the receiver is defined 
€=(4,/2,) sin (2nBp+8) (41) 


The angle of incidence of the ray trom the center or tne 
doublet image to the receiver intersecting the apparent 


pbettom is 
O=Z (11-1) Prbte (42) 


However, the upper and lower images in each doublet nave 
Sligntly different angles ot incidence with the apparent 


bottom 
0,= (241-1) Prdte, (45) 


The upper and lower angles at the recelver ror each aoubiet 


are 


€,=et (2,/2,) YCOS (2118 +8) (44) 
The individual angles of incidence are 
Oy-0- 7041+ (r,/r,) cos (2znptd) )+{(2,/r,) ycos (2B) (45) 


If the following variables are used in the coefficients as 


approximations ot the sine or the angles ot incidence 


ae 
a,71+— cos2nf 
r 
i ( 469 


P. 
b,-~— coc2nb 
Fo 


then eq. 38 becomes 


=e “Hane yb (47) 


R 


n, 


For each doublet image there iS an upper and lower image. 
The upper doublet 1s detined by A, With the upper image A. 
and the lower image A,. Likewise, the lower doublet b, has 


the upper image B. and the lower image BE.. 


A =&.e Hee nb,ya B= ee eee 
- n - 2 





A,= Re te ate Bane Ee oe a 
A, and AA detined by eq. 15 can be combined with eq. 438 
A,=%,cosh (nb,ya)e ™°* 8B =®.cosh(nb,ya) e72°* ue 
AA=8,sinh(nb,ya)e "=" = AB=&, sinh (nb,ya) e**" 
When combined, they yield 
£4 =tanh (nb_ya) 48 -+tanh(nb,ya) (50) 
Nn Qn 
The result is the full quadruplet equation 
N 25 
P=p,+}). £8 _coshb ei'*t *) e-i¢ 
DP, »y, I n (51) 


[e-2-JHsin (kr, ysin(2nB+8) ) -jtanhb cos (kr,ysin(2nf+d) ) 
-e**F4sin(kr,ysin(2nB-8) ) -jlanhb cos (kr, ysin(2nB-S) )] 


14 





Vile. RESULTS 

Three runs at dirtferent ranges were made to compare tne 
output or tne quadrupiet model with tnat or tne "URVEXT" 
program. For each run, the source angle is 1°, ana the weage 
angle is 3°. The ratio ot the sound speed in seawater to 
that of the bottom is 1.01695, and the ratio ot the density 
of seawater to that of the bottom is 0.7. This 18S equivalent 
to a Slow bottom consisting of a combination of clay and 
ooze; conditions very common 1n continental sheit regions. 
The runs are tabulated at 0.30° increments ot receiver angie 
from the surface downwards. The resultsS are snown 1n tabies 
1 through 3. At a recelver angle ot zero all three runs 
result in an amplitude of zero, as they snould. 

At the shortest range, Run 3, the resuits of tne 
quadruplet model differ signiricantly trom those or “URTEXKT . 
the percentage dirterence in amplitude 1S consistentiy apove 
10%. With the exception of the 3° receiver angie, the 
average percentage difference 1n amplituae tor the otner two 
runs 18S well below 10%. Tne tar tftield daoubiet approximation 
improves at longer ranges, and the approximations tor Ar 
values become more accurate. Also, the difference in phase 
angle is much greater at the shortest range; half the phase 
angles in Run 3 differ by greater than n/4, but in the other 
two runs, only the phase angle at a 3° receiver angle ditters 
by more than that amount. 


io 


Each of these data sets 1S accompanied by sampie piots 
of the quadruplet pressure amplitude and tne rerilection 
coetficient for each quadruplet. In each case, the plots 
were taken with a receiver angle of 0.90°. The magnitude ot 
the complex quadruplet pressure drops oft very quickly with 
the quadruplet number. It 1s negligible by tne 5th 
quadruplet. This 1S a reSuit otf the effect or the rapia 
(exponential) decay of the Raylielgh retlection coetticient 
with increased angle of incidence ot tne higner number or 
quadruplet. 

In practical terms, this means that the quadrupiet 
expansion becomes less accurate witn larger receiver and 
wedge angles. AS thneSe angieS I1ncreaSe, the ditterence in 
the angles of incidence pnetween the upper doublet sound patn 
and the jiower doublet sound patn wlil increase 
proportionately. Wnen these ditterences are ralsea to the 
exponential! power in the reriection coetricients, tne 
differences 1n amplitude between upper ana lower douplets at 
the receiver will be quite large. However, the tundamental 
assumption in the quadruplet expansSion 1s that the upper and 
lower doublet images are equivalent in amplitude, but 
opposite in phase. Therefore, the quadrupiet expansion 15s 


only vaiid tor very sma!il angies and iong ranges. 
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BETA=3,GAMMA=1,R1=1,R2=100,DENSITY RATIO=. 
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TABLE 1 


9,SPEED RATIO=1.0ne2 
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0.000092 | 8.4 iase | -23.9 i ae ' 
0.000175 ! 2.8 118.8 | -23.3 h42el } 
0.000243 : 1.25 419.2 on: te , 
0.000288 | 10.0 18.3 | -20.8 hog ; 
0.000306 | “ras eepi7uS | <16a6% een 
u.000295 ; 1.7 |16.2 j|-15.4 } 
0.000255 , 12.0 / 16.0 Hes ee f 
0.000190 | 5.0 ;25.4 log OP 
0.000112 } 13.9 ; 18.2 1.40 
0. COGi@ES | 295.0 oo 8s. . py OL. 6 , 


18 





AINEKATY FATANApHYNA Ja Apnzyrvuhoy 


Quadruplet Number 


FUBTOLFJa0Q9 uorzoaT Joy ysrapAey 


"@6° @="} 19d 


‘TH 





uny 


PUSTIII FZa0Q 


TABLE 2 
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DELTA | URTEXT 
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| | QUAD | As | URTEXD | QUAD | 
f AMP. AMP. AMP PHASE | PHASE |! PHasE | 
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APPENDIX A 


EEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEEES 
*S*QUADRUPLET EXPANSION 

*pquad.m 

*MATLAB 

*Quadruplet expansion ot tne acoustic 
pressure field in a wedge shaped ocean. 
*Michael Joyce 


SEEEEEEEEEEEEEEEEEEEEEEEEEEEESEE ES 


B=input('wedge angle in degrees='); 

G=input( ‘source angle from surtace in degrees='); 
D=1nput( receiver angle trom surtace 1n degrees= '}); 
rho=input('water to bottom density ratio='); 
cc=input('water to bottom sound speed ratio='); 
rl=input('range of source from apex='); 


r2=input('range ot receiver from apex='); 


N1l=£1x(90/B); 
B=B*pi/1380; 
G=G*p1/180; 


D=D*pi/180; 


SEEEEELEEEEEEEEEEEEEEEESEEEEEEEEES 


CALCULATE SCALING FACTOR 


GEEEEEEEEEEELEELEELEEEEEEEEEEEEEEEESD 


tb=tan(B); 
Lf ce<i1, 
tl=acos(cc); 
t2=sin(tl); 
else 
tl=acos(l/cc); 
t2=tan(tl); 
end 
kKI=2*t2*tb; 


k=pi/kl; 


SEKEBEEKBSBEBEESEREBESESBEEDEEESSESEDS 
*%CALCULATE CONSTANTS AND PRIMARY DOUBLE’ 


EEEKEBEEBREBEKREEEEEEEEEEEEEEEEEESESES 


al=2*(1/rho)/(sqrt(ce 2>i9oe 
r=sqrt(rl 2+r2 2-2*ri*r2*cosw ), 
r3=abs(r2-rl); 

r4=ElArZ7 eS 25 

mu=r2/abs(r2-rl1); 

q=k*r1*G*D; 

pl=(2/r)*sin(q); 


f=0; 
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EEEEEEEEEEEEEEEEEEEEEEESEBEEEEEES 
%QUADRUPLET SUMMATION 


SEEEEESESEEESEEEEESEEEESEEEEEEEEE 


ber n=1:1:Nl1, 
s(n)=n; 
th(n)=2*n*B; 
d({(n)=k*rl*D*¥sin(th(n)); 


prmeonm)=k*ril*mu*( l-cos(thi(n))); 


EEEEEEEESEESEEEEESEEESEES 
aVeorusellONn COEFFICIENTS 


SEBEEEBEEESESEEBEEESEEEEES 


Milnes err ler t+h2 2-241 i*r1 2*cos(thin)+D)); 
Been -Saqretrl 2+r2° 2-2*rl*r2*cos(th(n)=D) ): 
Both }=bi/ rutin) ; 

Bony Girl (nn); 

Geen =—l1rrotn jy cos(th(n) ); 

pin) =r6(n)*cos( thin) ) ; 

rr(n)=-B*ai*(n 2); 

R(n)=-exp(rr(n)); 

a(n)=n*al*D*al(n); 


b(n)=n*al%G*bl(n); 


SEEEEEEEESEEEEBEEBESEEEEES 
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%PRESSURE 


SESEEESEEBEREEEEESEEEEEES 


ul(n)=cosht bin); 
u2(n)=exp(-j)*phi(n)),; 


u3(n)=(-2/r)*R(n)*ul(n)*u2(n); 


EELEELECEEEEEEEEEEEEEEEES 
%UPPER IMAGE 


EESCEEBEEEEEEBEEEEETEESDEES 


V1I(n)=s2n( tne + Dy 
v2(n)=k* rl *GAviien ) > 
v3(n)=sin(v2(n)); 
v4(n)=cos(v2(n)); 
v5(n)=-a(n)-37*mur*d(n); 
vo(n)=exp(v5(n)); 
v7(n)=tanh(b(n)); 


ve(n)=avb6(n)*(v3(n)-J* v/(n)*v4a4(n)); 


ELEEEEEEEEELEEEEEEESESEEEE 
SLOWER IMAGE 


SECEEEEEESEEEEEEEEEEEEEES 


wl(n)=sin(th(n)-D); 


w2(n)=k*r1l*G*wl(n); 
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w3(n)=sin(w2z(n)); 
w4(n)=cos(wz(n)); 
w5(n)=a(n)+j3*mu*d(n) ; 
w6(n)=exp(w5(n)); 


w7(n)=w6(n)*(w3(n)-j*v7(n)*w4(n) ); 


SEEEEEEEEEEEEREEEEEESEEES 
%*SUMMATION 


KEEEKEERESEKRBEEBRESBEED EEE SESE 


p(n)=u3(n)*(v8(n)-w7(n)); 
P=t+pi(n); 
pr(n)=real(p(n)); 
pim(n)=1mag(p(n)); 
pz(n)=sqrt(pr(n) 2tpim(n) 2); 
Env = =i “RW Ty; 
end 

Peers, p2);gr1a; 

xlabel('Quadruplet Number'); 

ylabel( 'Magnatude ot Quadruplet Pressure'); 

title('Quadruplet Pressure'); 

pause; 

mlotts,rm)-grid: 

title( Rayleigh Reflection Coetricient ); 

xiabel( 'Quadruplet Number'); 


yiapel( Coerticient'); 


pause; 

ps=pl+t; 
ph=imag(ps)/reai(ps); 
phl=atan(ph); 
phase=180*phl/pi 

f=f; 

pl=pl; 

PS=ps; 

P=sqrt(real(ps) 2t+timag(ps) 2) 


end; 
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